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Abstract. Spin dependent single electron tunneling in a ferromagnetic double junction 
is investigated theoretically in the limit of incoherent sequential tunneling. The junction 
consists of a small nonmagnetic metallic grain with discrete energy levels, which is connected 
to two ferromagnetic electrodes. We have developed a new theoretical model taking into 
account charge as well as spin degrees of freedom. The model allows us to investigate new 
phenomena such as spin accumulation and spin fluctuations. 
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1 Introduction 

Single-electron tunneling in mesoscopic double-junctions has been extensively stud- 
ied in the last decade both experimentally and theoretically [Q. When the central 
electrode (grain) is small enough, then discrete charging of this electrode leads to 
Coulomb blockade of electric current below a certain threshold voltage and to charac- 
teristic 'Coulomb staircase' at higher voltages. The interplay of ferromagnetism and 
discrete charging was studied only very recently |^ . It was shown that spin accumu- 
lation in the grain, together with discrete charging, is responsible for oscillations in 
the tunnel magnetoresistance (TMR), i.e., in the change of junction resistance when 
the magnetic configuration is varied from antiparallel to parallel alignment. In this 
paper we study the influence of discrete energy levels |^ of the grain on single-electron 
tunneling phenomena in ferromagnetic junctions. The studies are performed within 
an approach, which is a generalization of the orthodox method [1] for the case with 
many electron channels including spins. 



2 Description of the model 

We consider spin dependent transport in a double ferromagnetic nanojunction in the 
limit of incoherent sequential tunneling, where the orthodox tunneling theory is appli- 
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cable. Accordingly, the resistances of the left, R2cri and right, i?icr, barriers for both 
spin directions (cr =t, i) are assumed to be much larger then the quantum resistance 
Rq = V2e^ 

In the stationary state the electric current / is given by 

I^^Y. Y.^rTl^^N^,N^)+rTl^{N^,N^)]P{N^,N^), (1) 

N^.N^ r=± 

where the probability P{N-^, N^) to find in the grain and N^^ excess electrons with 
spin cr =1 and a —I, respectively, is calculated from the relevant master equation [Q. 
In the present model we operate in the two-dimensional space of states (N-^ i Ni), in 
contrast to the spinless orthodox method [1] which was confined to the one-dimensional 
space. The spin dependent rates of electron tunneling to (r — +) and off (r = — ) the 
grain through the j-th junction = 1, 2) are given by H] 



1 -I- exp T 



\kBT 



-1 



(2) 



where Ep denotes the Fermi energy and the summation runs over discrete energy levels 
Ei„ of the grain, which are assumed to be equally separated with the level spacing 
/S.E. The formula (2) is obtained for fc^T « /S.E, when the distribution function of Na 
among the levels of the island differs appreciably from the Fermi-Dirac distribution 
function The voltage drop Vj{N^, N^) on the j-th junction is given by 

y,(iVT,iv,)^(-i)^y^ + ^, (3) 

where N — N-\^ + Ni \s the total number of excess electrons on the grain, Cj is 
the capacitance of the j-th junction, C denotes the total capacitance of the grain C = 
Ci + C'2, while y, e and T stand for the bias voltage, electron charge and temperature, 
respectively. Finally, Ec is the charging energy, Ec = e^/2C. We have also assumed 
that the energy relaxation time is much shorter than the time between successive 
tunneling events and shorter than the spin relaxation time. Magnetic configuration 
of the junction can be varied by an external magnetic field and the rotation from 
antiparallel to parallel alignment is accompanied by a change in the junction resistance 
(in our case Rii, Rii Ri^)- 



3 Results and discussion 

We have calculated numerically the bias voltage dependence of the current /|| and 
flowing through the junction with parallel and antiparallel alignment of the electrode 
magnetic moments, TMR— (/|| — /||)//||, spin accumulation < M >=< — N ^ > 
and spin fluctuations rms(Af) = y/<AP>^^<TT>^ . The results are presented in 
Fig.l and Fig. 2 for various temperatures. For T = 2.3 K each curve in Fig.l has 
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Fig. 1 Voltage dependence of the current in the antiparallel configuration (a), the tunnel 
magnetoresistance (b), the spin accumulation (c) and the spin fluctuations (d) for different 
temperatures. The parameters are: Rif = 2 MQ, Rii = 4 Mf2, J?2t = 200 Mfl, R21 = 100 
Mn, Ci = 4.3 aF, C2 = 9 aF, Sc/fcs = 69.9 K and AE/ks = 34.8 K. 

additional steps which result from opening additional electron channels corresponding 
to discrete energy levels. These steps are clearly seen only for T <^ AE/ks, while 
for higher T they are smeared out. Spin asymmetry in the tunneling rates leads 
to imbalance of incoming and outgoing currents, which causes a change in the local 
chemical potential and in the number of electrons of a particular spin orientation. 
This, in turn, gives rise to spin accumulation (Fig.lc). One can note that in the 
first stage the average < M > increases with increasing bias voltage, which is due 
to increasing number of active magnetic channels with increasing V (see also spin 
fluctuations presented in Fig. Id). This process continues until 1/ « 20 mV, when it is 
more favorable to open a new charge channel. This, in turn, results in a reduction of the 
spin accumulation and spin fluctuations. At low T the shape of the spin accumulation 
curve for even number of electrons accumulated in the grain is different from the 
shape corresponding to an odd number of electrons (this is due to a difi^erence in the 
(A''^, Ni) - space of states for the cases of even and odd numbers). The effect of 
spin accumulation has a direct influence on the character of TMR (see Fig. lb, where 
the small steps correspond to discrete energy levels while the large-scale oscillations 
correspond to single electron charging). 

The situation analysed above corresponds to a system with a large di£Fcrcnc;c in 
the junction resistances. In the case of a symmetrical double junction, on the other 
hand, the /-F steps are not seen in the parallel configuration as shown in Fig. 2a 
(apart from small steps scon at low temperatures, which result from the discreteness 
of electronic structure). However, the value of TMR in the ohmic limit (for high T) 
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Voltage [mV] Voltage [mV] 

Fig. 2 Voltage characteristics of the current (a) and TMR (b) for the device with symmetric 
tunnehng resistances = i?2T = 102 Mfi, R\y = R21 — 51 Mfl in the parallel configuration. 
The other parameters are the same as in Fig.l. 



can be larger (compare Figs lb and 2b). In the antiparallel configuration, on the other 
hand, the junction ceases to be symmetrical, which leads to the oscillations in TMR 
with increasing V, seen in Fig. 2b. However, amplitude of these oscillations is quickly 
damped. 



4 Summary 

We presented in this paper calculations of electric current in ferromagnetic double 
nanojunctions with discrete electronic structure of the central electrode. The discrete- 
ness leads to additional structure in the /- V curves and in the voltage dependence 
of spin accumulation, spin fluctuations and TMR. This fine structure is clearly seen 
only for small T, while for higher T it is smeared out. When the spin-flip relaxation 
time is much longer than the time between successive tunneling events, the effect of 
spin accumulation leads to TMR which oscillates with increasing V (clearly seen for 
asymmetrical junctions). 

The paper is supported by KBN under Grant No. 2 P03B 075 14. 
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